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Abstract 

We study the equilibrium, oscillatory and transport properties of a chain of charged particles 
which interact with each other via the Coulomb and power-like repulsive interactions. Exact, 
analytical expressions for the energy of the ground state, interparticle distances, and vibration 
spectrum are derived, and the stability criterion for a chain with equidistant particles is established. 
We also consider the phenomenon of deterministic transport of a chain in a ratchet potential subject 
to a longitudinal alternating electric field. 
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I. INTRODUCTION 



A chain of coupled oscillators is a universal mathematical model that describes a large 
variety of phenomena in physics, chemistry biology and other sciences. Examples include 
energy localization and propagation, 1-3 discrete solitons, 4 ' 5 chemical reactions, 6 nonequilibri- 
um phase transitions, 7 ' 8 stochastic resonance, 9 ' 10 synchronization, 11 brain dynamics 12 and 
many others. A special case of that model, a chain of charged particles that interact via the 
Coulomb interaction, also has a variety of applications. In particular, if the charges have the 
same sign and are confined by a trapping potential, it describes the so-called one-dimensional 
Coulomb crystals, i.e., a crystallized structure of ions in a linear trap formed by external 
focusing fields. 13-15 If neighboring charges have opposite signs and a repulsive interaction 
exists between particles, then a such chain represents a relevant model for one-dimensional 
ionic crystals. It seems that those models can be useful also for the study of one-dimensional 
crystals within carbon nanotubes. 16,17 

During the last decade, much attention was devoted to the transport properties of par- 
ticles in thermal ratchets, i.e., spatially periodic stochastic models exhibiting noise- induced 
macroscopic transport. 18-27 The ratchet effect or stochastic transport in asymmetric poten- 
tials is important for biological motors, 28 ' 29 separation techniques, 30 surface smoothing 31 and 
vortex density reduction, 32 to name a few. Most studies of ratchets, including deterministic 
ones, 33 ' 34 i.e., ratchets without internal and external noises, have been performed for a single 
particle. However, examples of coupled particles 35,36 and vortices 37 show that their collective 
interactions can significantly modify the transport properties from the single particle case. 
To our knowledge, the transport properties of a chain of strongly coupled charged particles 
in deterministic ratchets have not been studied at all. 

In this paper, we study the equilibrium, oscillatory and transport properties of a chain 
where the neighboring particles have the same mass and opposite charges. Particles interact 
with each other via both the Coulomb interaction and a repulsive interaction that depends 
on the interparticle distance I r (xi = Xi(t) is the coordinate of the ith particle, 

r > 1). In Sec. II, we consider the ground state of such a chain depending on the exponent 
r. Three different states of a chain, namely, a periodic chain with equidistant particles, a 
periodic chain with non-equidistant particles, and a chain with infinite period, are analyzed, 
and the interparticle distances and the chain energies as functions of r are calculated. In 
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contrast to intuitive expectations, we show that there is a critical value of the exponent, 
r = r 2 ~ 2.799, such that for r > r 2 the minimum energy occurs for a chain with equidistant 
particles, and for r < r 2 the minimum energy occurs for a chain with infinite period. 

Propagation of longitudinal waves in a chain with equidistant particles, its stability with 
respect to small displacement of particles from their equilibrium positions, and forced vibra- 
tions of particles under action of an alternating electric field are studied in Sec. III. We derive 
the set of equations that describe small oscillations of the particles, find the phonon spec- 
trum and its asymptotics near the center and boundary of the Brillouin zone, and calculate 
the amplitude of forced vibrations. We show that the spectrum consists of one branch; on 
the boundary of the Brillouin zone it has a non-analyticity which is caused by the long-range 
action of the Coulomb interaction, and a chain is stable only if r > r\ ~ 2.779. 

The deterministic transport of a chain is studied in Sec. IV. We consider the case when the 
ratchet potential is charge independent, and a chain is driven by a longitudinal alternating 
electric field. We find the average velocity of a chain and analyze in detail its dependence on 
the parameters characterizing the interparticle interactions, ratchet potential and alternating 
electric field. Concluding remarks are contained in Sec. V, and some technical details are 
given in the Appendixes. 



II. GROUND STATE 



We consider a chain of particles which have the same mass M and any neighboring 
particles have opposite charges, i.e., q i+ i = —qi and \qi\ = q (the index i labels the particles). 
We assume that in addition to the Coulomb interaction a repulsive interaction exists between 
the particles. In analogy to ionic crystals, 38 we choose a repulsive interaction that depends 
on a power of the interparticle distance \xi — Xj\, and write the total interaction energy of 
a chain in the form 

^sVV^fV (2 ' 1) 

Here b is a dimensioned parameter characterizing the strength of the repulsive interaction, 
and the primes on the summation signs indicate the absence of self-interaction terms, i.e., 
i 7^ j. To prevent the unbounded contraction of a chain, we require r > 1. 

If a chain is periodic in equilibrium (see Fig. 1), then the coordinates of the positive and 
negative charges are written as np + x\ and np + a + x~, respectively. Here p is the chain 

3 



period, n (n = 0, ±1, ±00) numbers the chain cells which contain two particles, a is the 
interparticle equilibrium distance in the cells (a < p), and x+ = x+(t) and x~ = x~(t) are 
the displacements of positive and negative charges from their equilibrium positions. With 
these definitions Eq. (2.1) reads 
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1 
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(n-m)p + x+ -x+\ r 
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\(n-m)p + x n -x; 



(2.2) 



(n — m)p — a + rc+ — x m \ r 

where summation on n and m is taken over all integers. Note also that the particles cannot 
move past each other, and therefore the conditions a > x+ — x~ and p — a > x~ — x^ +1 must 
hold for all n and t. 

In equilibrium x+ = x~ = 0, and Eq. (2.2) is reduced to the form 
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(2.3) 



where N is the number of particles in a chain (N — > 00), ct = a/p (0 < ct < 1), and 
C( r ) — D^Li n ~ r is the Riemann zeta function. Taking into account that the energy W^ t 
is invariant with respect to the replacement of a by 1 — a, we conclude that such chains 
are equivalent. Therefore, the permissible values of a can be restricted by the condition 
< a < 1/2. 

Defining the dimensionless parameter (3 — p/l (0 < (3 < 00) and the length scale / = 
( 6 / 5 2)i/(r-i) > we can rewrite E q. (2.3) as wj$ = N(q 2 /l)U (a, (3), where 



U (a,P) = 



it 



+ 



C(0 



2a/3 /3 ^ n(n 2 - a 2 ) (3' 



+ —Y^ - (2 A) 

2(3 r z ^\n + a\ r K ' J 

is the dimensionless chain energy per particle (for another representation of Uo(a,(3) see 
Appendix A). The equilibrium values of a and (3 are defined by the set of equations 
dU (a, f3) /da = and dU (a, f3)/df3 = 0, which are written as 



1 2a ^ n r n + a 

2a~ 2 f3 ~ Thi (n 2 - a 2 ) 2 ~ W r \ \n + a\ r + 2 = ' 



1 a^ ~ 1 r_ ( 

2a(3 2 + (3 2 ^ n(n 2 - a 2 ) ^ Q{V) 



(2.5) 



+ (2 - 6) 

Analysis shows that these equations have three solutions: a = a s (r) and (3 = (3 s {r) 
(s = 1,2,3). Below we find a s (r) and (3 s (r), and calculate the chain equilibrium ener- 
gies U s (r) = U (a s (r), f3 s (r)). For illustrative purposes we introduce also the parameters 
7 s (r) = a s (r) (3 s (r) , which define the equilibrium values of the parameter 7 = af3 = a/l. 



A. First solution 

For a — 1/2 (p = 2a), i.e., for a chain with equidistant particles, the series in Eqs. (2.4) - 
(2.6) are either standard series or can be reduced to them. Indeed, using standard series, 39 
we find with linear accuracy in e = 1/2 — a (e — > 0): 

^w^f^ 2 - mi3) - 6]£ - (27) 

E 1 " + L = 2 ( r + 1 K 2 " +2 - l)C(r + 2) £ , (2.8) 
00 j 

EjjTOj =4(21n2-l)-8(3-41„2) £ , (2.9) 

E^p-^-DCM. (2.10) 
According to these formulas, [/ (l/2,/3) is written as 

C/o(lA/3) = -|ln2+^ C(r), (2.11) 
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at a — 1/2 Eq. (2.5) is satisfied identically, i.e., a±(r) = 1/2, and Eq. (2.6), which at a = 1/2 
is equivalent to the equation dU (l/2, f3)/d(3 = 0, yields 

ftW = 2 ( I S 1 ) 1/< " ) ' (2 ' 12) 

Substituting f3 = /3i(r) into Eq. (2.11), we find for the equilibrium chain 

r- 1 / In 2 \ 1/(r_1) 

w= - w tQ) ' (2 - 13) 

[Note that U\{r) is the minimum value of Uq{1/2,(3) only if r > 1.] 

The plots of the functions /3i(r), 7i(r), and £/i(r) are shown in Figs. 2, 3 and 4, respec- 
tively. To gain more insight into the behavior of these functions near r = oo and r = 1, we 
define the limiting values of those functions as r — > oo and their asymptotics as r — > 1. Tak- 
ing into account that lim x _> 2; x = 1, lim x _ ( .o(l + 2;) 1 ^ x = e > C(°°) = 1> an d C( r ) ~ (r — 1) _1 as 
r -> 1, we obtain from Eq. (2.12) /3i(oo) = 2 and/3i(r) = 2e(eln2) -1 / e -> oo as e = r-1 -> 0, 
and from Eq. (2.13) E/i(oo) = - In 2 and U x {r) ~ -e _1 (e In 2) 1+1 / e -> as e -> 0. Thus, if 
r — > 1 then the functions /3i(r) and £/i(r) go extremely fast to infinity and zero, respectively. 
From a physical point of view, this behavior of /3i(r) and Ui(r) is caused by the unlimited 
growth of the total repulsive energy of a chain as r — > 1. 

B. Second solution 

If (3 7^ oo, then we eliminate /3 from Eq. (2.5) with the help of Eq. (2.6) and obtain the 
equation 



r 



a 2 ~l (n 2 — a 2 ) 2 J \ ^\n + a 

( 1 9 2 1 \^n + a 

~U + a hi < n2 - « 2 ) / n \n + a\ r+2 = ' 

(2.14) 

which defines the equilibrium values of a. According to the results of the previous subsection, 
Eq. (2.14) always has the solution a = 1/2, i.e., a\{r) = 1/2. To determine if Eq. (2.14) has 
solutions on the interval (0, 1/2), we find first the asymptotics of its left-hand side, denoted 
by L r (a), as a — > and a — > 1/2. If a — > 0, then 

L r (a) ~ [a- 2 -4C(3)a][a- r + 2C(r)] 

-[«- 1 + 2C(3)a>- (r+1) , (2.15) 
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and we obtain different asymptotics of L r (a) for different values of r: 

' -6C(3)a-( r - 1 ), r > 3, 
L r (a) ~ { -2((3)a- 2 , r = 3, 
+4C(r)a" 2 , r < 3. 

Using Eqs. (2.7)-(2.10), as a -> 1/2 we obtain 

r+3. 



(2.16) 



where 



L r (a) ~ -2 r+d C(r)F(r)e, 



F(r) = 4(1 - 2— 2 )(r + 1) l^^y^ - 7C(3) 



(2.17) 



(2.18) 



is a monotone increasing function of r, and F(r) = at r = 77 2.779. Therefore, in the 
left vicinity of the point a = 1/2, the conditions L r (a) > and L r (a) < hold for r < r 1 
and r > ri, respectively. Using this fact and Eq. (2.16), we conclude that if 77 < r < 3, then 
L r (a) near the points a = and a — 1/2 has different signs, and so Eq. (2.14) has at last 
one solution a 2 (r) G (0, 1/2) satisfying the conditions ai2(r) - >■ 1/2 as r — * ri and a2(r) — > 
as r — > 3. Numerical analysis shows that such a solution is unique, and Eq. (2.6) yields 



2CW + E 



1 



n + a 2 (r) 



+ 2« 2 2 (r)E 



i 

r-l 



(2.19) 



a 2 (r) ' '"'^ vy ^ 1 n[n 2 -ai(r)] 

Since a chain is transformed to a chain with equidistant particles as r — > ri, the equalities 
/M r i) = A( r i)5 72( r i) = 7i( r i), and U 2 (r 1 ) = Ui{ri) hold. If r — > 3, then according to 
Eq. (2.19) p 2 (r) ~ V3/a 2 (r) -»• oo, and we obtain 7 2 (3) = v 7 ^, ^(3) = -1/Zy/Z. The plots 
of the functions (3 2 (r), 72 (r), and U 2 (r), which are defined for r e [ri,3], are also shown in 
Figs. 2-4. As is seen from Fig. 4, U 2 (r) > U\(r) for 77 < r < 3, i.e., the equilibrium chain 
with equidistant particles is energy-wise preferred compared to the equilibrium chain with 
non-equidistant particles. 



C. Third solution 

A further solution of Eqs. (2.5) and (2.6) corresponds to the case when the period of the 
chain is infinite and the distance between particles in the cells is finite, i.e., when a — > 0, 
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(3 — > oo, and a(3 = 7. In this case 

Uo (a,P) = U ( 1 ) = -±- + 1 ^, (2.20) 

Eq. (2.6) is satisfied identically, and Eq. (2.5) is reduced to the equation dUoij)/d^ = 0. Its 
solution 7 = 73(7"), 73(7") = r 1 / (r ~ 1 ), defines the equilibrium distance between the particles, 
and the equilibrium energy per particle is given by 

r - 1 /lxVtr-i) 
W = -— Z • (2-21) 



2r 

The functions 73 (r) and U%{r) monotonically decrease as r grows from 1 to 00 [C/ 3 ~ (r — 
l)/2e -> as r -> 1, C/ 3 (3) = -l/Zy/Z, U 3 (oo) = -1/2, 73 (1) = e, 73 (3) = y/Z, 73(00) = 1], 
and their plots are shown in Figs. 3 and 4. Comparing the functions £/i(r), U 2 (r), and t/ 3 (r), 
we conclude that for r > r 2 [r = r 2 ~ 2.799 is the solution of the equation C/i(r) = ^(r)] 
the minimum energy occurs for a chain with equidistant particles [p = 2a, a = lji(r)], and 
for r < r 2 the minimum energy occurs for a chain with infinite period [p = 00, a = £73(7")]. 
Note also that a chain with non-equidistant particles has a slightly larger energy than other 
chains. 

Next we study the vibration properties of a chain with equidistant particles. 

III. VIBRATION PROPERTIES 
A. Dispersion relation 

Assuming that the condition \x^ \ <C a holds for all particles, we expand the chain energy 
Wint in a power series of x^/a with quadratic accuracy: Wi n t = W^j + W^t + ^inh where 
the upper indices indicate the power of Straightforward calculations based on Eq. (2.2) 
with p = 2a and arbitrary value of a yield wj$ = N(q 2 /l) U (l/2, 2a//), wj$ = 0, and 



W int = 7E' 5 2(«-m) \{x+ - X^f + (x n ~ X 



ml 



where 



+ 1 -Y, B ^n-m)-l{xt-X- m )\ (3.1) 



(-l)"2g 2 6r(r + l) 
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For the case considered here that a has an equilibrium value, i.e., a = lji(r), we represent 
B n in the form B n = Mtt 2 B n , where fi 2 = q 2 /Ml 3 , 



B 



n 3/ \ 



o(-l) w , (r + l)ln2 

1 77- 1 3 C( r )l n | r+2 



(3.3) 



and according to Eq. (2.12) 7l (r) = [r((r)/ ln^ 1 ^" 1 ). 

From Eq. (3.1), we obtain the system of equations, Mx^ + d W^i /dx^ = 0, that describe 
free longitudinal vibrations of particles near their equilibrium positions in the harmonic 
approximation, in the form 

X n + Q2 2 ^2(n-m)-l(^n ~ ^m) 

+ n 2 ]T 5 2(n _ m) (x+ - x+) = o, 

(3.4) 

+ ^ Z! ^2(n-m)+l(^ n ~ &m) 
m 

+ fi 2 ]T B 2 (n-m)(x~ - X m) = °- 

Since n ranges over the integers, this system contains infinitely many coupled equations. We 
seek wave solutions, which in the complex representation have the form 

x n = x ± Gxp(i2nak — iut). (3.5) 

Here x + and x_ are the complex amplitudes, k is the wave number belonging to the Brillouin 
zone, i.e., to the region [— 7r/2a, n/2a] of reciprocal space, and uj is the wave frequency. Sub- 
stituting solution (3.5) into Eqs. (3.4) and using the definition k = ak/n (k e [—1/2, 1/2]), 

oo 

D = 2Y,B P , (3.6) 
P =i 

oo 

P(k) = 2^5 2p COs(2j97T/€), (3.7) 

p=l 

= X/ ^2p-i exp(— i2piiK), (3.8) 
p 

we obtain the system of equations for x + and x_ 

x + [-^ 2 + D - P(k)] - x_S(re) = 0, 

(3.9) 

-x + S*(re) + x_[-£ 2 + £> - P(«)] = 0, 
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where uj = oj/VL is the dimensionless wave frequency, and the asterisk denotes complex 
conjugation. This system has non-trivial solutions, i.e, propagation of phonons along a 
chain can occur, if its determinant equals zero, 41 

[u; 2 -D + P(k)} 2 -\S(k)\ 2 = 0. (3.10) 

Solving Eq. (3.10) with respect to u, we find the dispersion relation for phonons: Ui^k) = 
Qui^k), where 

0)l )2 («) = y/D-P(K)±\S(K)\, (3.11) 

and the plus and minus sign in front of ^(/c)! correspond to Cj\{k) and u) 2 («), respectively. 



B. Analysis of the spectrum 



Since u)i )2 (— k) = cj lj2 (/t), we can restrict the permissible values of the dimensionless 
wave number k by the condition k G [0, 1/2]. According to Eq. (3.11), in this interval the 
dimensionless phonon spectrum consists of two branches, Cj\{k) and u) 2 (k;). It is easy to see 
that on the interval [0, 1] the spectrum is given by 



Cj{k) 



&i(k), 0<k<1/2, 
tD 2 (l-K), 1/2<K<1, 



(3.12) 



and u>(k) is continuous at k — 1/2. Indeed, for k = 1/2 Eq. (3.8) yields 5(1/2) = 
^l p {— l) p -B 2p -i. Replacing in this formula p by — p + 1 and using the condition B n = B_ n , 
we obtain S{l/2) = - Ep(-l) p #2 P -i- This implies that 5(1/2) = 0, and Eq. (3.11) yields 
d>i(l/2) = d) 2 (l/2). Hence, the phonon spectrum on the interval [0,1] is represented by 
the single branch uj(k). Note that the region [—1,1] is the dimensionless Brillouin zone 
corresponding to a chain with a primitive cell which contains one particle. 

To study the behavior of the phonon spectrum £>(«), we use Eq. (3.3) to represent 
Eqs. (3.6)-(3.8) in the form 

1 



D = 
P{k) = 
S(k) = 



7 3( r ) 
1 

1M 



2(r + l)ln2 C( ^. + . 2) - 3((3) 



(r + l)ln2 

(r + l)ln2 

2C(r) 



C(0 

$ r+2 (/t) + $ 3 (k) 
^ r+2 (/t) - * 3 (k) 



(3.13) 
(3.14) 
(3.15) 
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where 



p=l p 



OO 



cos|p-^_ 
v; p t! (2p-l) s 

[Recall that the functions P(k) and S(n) are defined for |/c| < 1/2.] Unfortunately, the 
analytical summation of the series in Eqs. (3.14) and (3.15) can be performed only for 
special cases, such as k = or 1/2. We calculate u>(k) only in the centre (k = 0) and 
on the boundary (k = 1) of the Brillouin zone which, according to Eq. (3.12), requires us 
to determine P(0) and 5(0). Using the definition of the Riemann zeta function and the 
standard series Y%Li{2p — l)~ s = (1 — 2~ s )((s), we obtain 

P(0) = -4rr 
1 ' 2 7l 3 (r) 

and D - P(0) = S(0) = F(r)/2 7l 3 (r), where F(r) is defined by Eq. (2.18). If F(r) > 0, i.e., 
r > ri, then 



2-(r + l)ln2^±^ + C(3) 



(3.18) 



u>(0) = 0, ^(l) = v /F(r)/ 7l 3 (r). (3.19) 

Note that Eqs. (3.9) yield x + = x_ for Co = 0, and x + = — x_ for Co = Co(l). This means that 
neighboring particles in a chain perform acoustical vibrations at k — and optical ones at 
k = 1. 

To find for r > r\ the asymptotics of a)(/c) as k — > and k — > 1, we need to calculate the 
series (3.16) and (3.17) with quadratic accuracy in k as k — > 0. This is not a problem for 
the series $ r +2(ft) and \E , r+2 (fi;) that arise from the repulsive interaction of particles, because 
expanding the cosines in a power series of k to an accuracy of k 2 we immediately obtain 

*r+2(«) = C(r + 2)-27r 2 C(r) K 2 , (3.20) 
* r+2 («) = (1-2— 2 )C(r + 2)-y(l-2-)C(r)^ 2 . 

(3.21) 

As to the series <&z{k) and that arise from the Coulomb interaction, their calculation 

by this method is impossible since divergencies occur. To avoid this difficulty, we use their 
integral representations obtained in Appendix B. These representations have the advantage 
that the double integrals in Eqs. (B3) and (B6) can be easily calculated with arbitrary 
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accuracy as k — > 0. Specifically, if we use the approximate formulas In [2 sin(y/2)] my and 
ln[tan(y/2)] ps ln(y/2) (y < 1), then Eqs. (B3) and (B6) yield 

$ 3 («) = C(3)-7r 2 [3-21n(27r/t)]/t 2 , (3.22) 
*s(«) = ^C(3)-y (3-21n^) K 2 . (3.23) 

Substituting series (3.20)-(3.23) into Eqs. (3.14) and (3.15), we find the main terms of 
the asymptotic expansion of cu(k) as k — > 0, 



(!)(/«) ~ 7T, 

and as k — > 1 , 



(r-l)ln2 

\ 1 3 A «, 3.24 



^ ~ + a>(lh?(r) (1 " K)2{ ln(1 " (3 - 25) 
According to these asymptotics, the long- wave part of the vibrational spectrum (for k« 1) 

has the usual form, i.e., Cj{k) oc k. On the boundary of the Brillouin zone, the function u(k) 

has a local minimum, cD'(l) = duj(K) / dn\ K= i = 0, and the short-wave part of the spectrum 

(for k ~ 1) exhibits a non-analytic dependence on k. As Eqs. (3.22) and (3.23) show, 

the non-analyticity of lv(k) at k — > 1 is caused by the long-range action of the Coulomb 

interaction. 

The dependence of Co(k) on k for different values of r is shown in Fig. 5. As is seen 
from this figure, cD(/c) is a non-monotonic function of k which has a unique maximum on 
the interval (1/2,1). As r decreases, its local minimum d>(l) also decreases, and u>(l) > 
for r > 7*1 . Since the inequality cD(l) > implies that cD(/c) > for all k ^ 0, the condition 
r > ri is the stability criterion of a chain with respect to small displacements of particles 
from their equilibrium positions. Note that though for r\ < r < r 2 a chain with infinite 
period is energy-wise preferred compared to a chain with equidistant particles, it is the latter 
which is stable. If r < ri, then a region exists on the interval [0, 1] where Co 2 (k) < 0, and 
therefore such chains are unstable. Finally, if r = r±, then cD(l) = 0, cD'(l) = — oo, and the 
asymptotic of u(k) as k — > 1 has the form 



~ ^ry (! " K )Vl ln ( x - ( 3 - 26 ) 

In this case the optical displacements of particles do not lead to an appearance of restoring 
forces. Therefore, to determine if a chain is stable or not at r = r 1; it is necessary to go 
beyond the harmonic approximation. 
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C. Forced vibrations 



We briefly consider also the forced vibrations of particles under the action of the longi- 
tudinal alternating electric field E(t) = Ecosu> t, where E and u are the amplitude and 
frequency of the alternating field, respectively. In such field the electric force qE cos uj t acts 
on positive charges, and the force — qE cos uiot acts on negative charges. If the friction force 
— Xx^ (A is a damping coefficient) also acts on each particle, then the forced vibrations of 
particles are described by the following system of motion equations: 

a£ + 2Q x x+ + tt 2 J2 B 2 ( n -m)-i(x+ - x~) 



(3.27) 



+ tt 2 B 2 ( n -m)(Xn ~ ^m) = AcOSU t, 

(n-m)+l{x n X m ) 
m 

+ fi2 Yj B 2(n-m)(x~ ~ X~) = A COS Ulftt 

(Q\ = A/2M, A = qE I'M). Symmetry of these equations suggests that we seek their solution 
in the form x^ = x(t) and x~ = —x(t). Taking into account that Co 2 {l) = 2J2^L\ B 2p -i, in 
this case Eqs. (3.27) are reduced to the ordinary differential equation for a forced harmonic 
oscillator 

x(t) + 2Q x x(t) + uj 2 {l)x{t) = A cos w t. (3.28) 

Its steady-state solution, which describes the forced vibrations of particles, has the form 
x(t) = xo cos(a;oi — <f), where 

x = , A (3.29) 

v /[^(l)-c 2] 2 + 4^I 

is the amplitude displacement of particles from their equilibrium positions, and tan</? = 
2f2 A c<;/[cJ 2 (l) — lu 2 ]. According to Eq. (3.29), if > \/2Q\, then the dependence of x 
on the driving frequency ooq has resonance character; ooq = [uo 2 (\) — 2Vt 2 x } 1 ^ 2 is the resonance 
frequency, and x = A/2VL\[uo 2 (l) — Q 2 ] 1 / 2 is the resonance amplitude. 



13 



IV. DETERMINISTIC TRANSPORT 



In this section we study the transport properties of a chain whose particles are driven by 
the longitudinal alternating electric field E(t) = Eh(t), 



h(t) 



1, 2kr <t< 2kr + r, 

(4.1) 

-1, 2kr + r < t < 2kr + 2r, 



and interact with the sawtooth (ratchet) potential V(x) that induces the force field f(x) 
—dV(x)/dx = fog(x), where 



g{x) = < 



cl (d — c), md + c < x < md + d, 

(4.2) 

— 1, md < x < md + c. 



Here k and m run over all integers, r is the half-period of E(t), fo = \ min/(a;)|, and d is 
the period of V(x) (see Fig. 6). The dynamics of particles in such a chain is described by 
the system of equations (3.27) in which the right-hand sides Acosuot and — Acosu t must 
be replaced by Ah(t) + Rg(x+) and —Ah(t) +Rg(x~) (R = f /M), respectively. To simplify 
the problem, we assume that the condition a = Ld (L is a natural number) holds, so that 
the equilibrium positions of all particles coincide with the minima of V(x). In this case, 
under the action of the electric field E(t) all positively charged particles and all negatively 
charged particles move identically, i.e., x+ = w(t) and x~ = u(t), and we obtain 

w + 2£l x w + oo 2 (l)(w - u) = Ah(t) + Rg(w), 

(4.3) 

u + 2Q x u + u 2 (l)(u -w) = -Ah(t) + Rg{u). 

Due to the periodicity of h(t) and g(x), the steady-state solution of Eqs. (4.3) satisfies 
the conditions 

w(t + r) = u(t) + Kd, u(t + r) =w(t) + Kd, (4.4) 
where K is a whole number. It follows from these conditions that 

w(t + 2t) - w(t) = 2Kd, u(t + 2t) - u(t) = 2Kd, (4.5) 

i.e., for the period of h(t) all particles are displaced the same distance A = 2Kd. This fact 
permits us to define the average velocity of a chain as 

v ^ ( 2r)- m(0)=K 4 
2t t 
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To find v, we need to calculate A and in the general case solve Eqs. (4.3). These equations 
have the important feature that we need to consider only the ratchet potentials with c > d/2. 
Indeed, the transformation c — > d — c leads to the inverted potential V in (x) = V(— x) which 
gives rise to the reduced force field gi n (x) = —g(—x). According to Eqs. (4.3), in such a 
potential the displacements w in (t) and u in (t) of positive and negative charges from their 
equilibrium positions are governed by the equations of motion 

w in + 2Q x w in + Lu 2 (l)(w in - u in ) = Ah(t) - Rg(-w in ), 

(4.7) 

U in + 2Q X U in + UJ 2 (l)(u in - Wi n ) = -Ah(t) - Rg(-Ui n ). 

Comparing these equations with Eq. (4.3), we obtain w in (t) = —u(t) and u in {t) = —w(t). 
This means that if for c > d/2 the solution of equations (4.3) is known, then their solution 
for c < d/2 is known too. Using Eqs. (4.5) and (4.6), we find the chain velocity in the 
inverted potential v in = —v, i.e., in the potentials V in (x) and V(x) a chain drifts in the 
opposite directions with the same velocity. If c = d/2 then V in (x) = V(x), both conditions 
Vi n — —v and = v hold, and Vi n — v — 0. Hence, to study the transport properties of a 
chain, we indeed can consider only those ratchet potentials for which c > d/2. In this case, 
a chain can drift only along the positive direction of the x axis and so K > 0. 

Next we study the deterministic transport of a chain in the overdamped limit where 
Eqs. (4.3) are reduced to 

2VL x w + uj 2 {1){w -u)= Ah(t) + Rg(w), 

(4.8) 

2fl x u + u 2 (l)(u -w) = -Ah(t) + Rg(u). 

Since our aim is to find the chain velocity v, we need to solve Eqs. (4.8) only for t G (0, 2r). 
To simplify the problem further, we assume that the conditions L » i( » 1 and 1 < A/R < 
c/ (d — c) hold. The former implies that the interparticle distance a is much larger than the 
chain displacement A, which, in turn, is much larger than the potential period d. The latter 
shows that for h(t) > (h(t) < 0) positive (negative) charges can overcome the potential 
barriers of V(x) and they move along the positive direction of the x axis, whereas negative 
(positive) charges stay in the minima of V(x). After h(t) changes sign, the particles at 
rest begin to move, and the moving particles move to the appropriate minima of V(x) and 
stop there. Since the maximal distance that the moving particles travel after h(t) changes 
sign does not exceed c, i.e., this distance is much less than A, we assume that the moving 
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particles reach the corresponding minima of V(x) immediately as h(t) changes sign. Within 
this approximation, the dynamics of positive and negative charges is described separately, 
and Eqs. (4.8) yield u(t) = and 



2fl x w + uo 2 (l)w = A + Rg(w) 
HO) = -A/2, w(r) = A/2] for t e (0,r), and w(t) = A/2 and 



2tt x u + u 2 (l)(u - A/2) = A + Rg(u) 



(4.9) 



(4.10) 



[u(r) = 0, u(2r) = A] for t e (r, 2r). 

To calculate A we can use either Eq. (4.9) or Eq. (4.10). We take Eq. (4.9) as our starting 
point and write it in the form 



w + 5w = < 



02, md + c < w < md + d, 
9i, md < w < md + c, 



(4.11) 



where 5 = u 2 {l)/2Q x , 0i — (A — R)/2Q X , 9 2 = [A + Rc/{d - c)}/2Q x , and m = -K, -K + 
1,...,K — 1. The solution of Eq. (4.11) depends on A and actually defines A, because 
u>(0) = —A/2 and w(t) = A/2. We use this feature to derive the following equation (see 
Appendix C) 



-Mil 



(20^ \ r m-A*_c + ; 

I 2d5 I 2d5 d , 



\ 2d5 J \ 2dS d , 

( 29 2 + A5 _ c \ ( 29 2 - A5 \ 

{ 2d5 d j [ 2d5 j 

x- 



, , 29 2 -AS c \ 29 2 + A5 ' 
[ 2d5 d [ 2d5 



(4.12) 



for the instant t = t 2 x that corresponds to a displacement of positive charges of A = 2Kd. 
An increase of r leads to the step- wise increase of t 2 K and K, and from Eq. (4.12) we obtain 



max K = min 
where [x] denotes the integer part of x. 



9 x + {d- c)S 
d5 



02 

d5 



(4.13) 
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Replacing in Eq. (4.12) t 2 K by r, we obtain the approximate equation that defines A as 
a function of r. For 5t 3> 1 it yields K maxi^, and so in this case the chain velocity v is 
given by 

2d 

v = — max If. (4.14) 

T 

It follows from this formula that v — > as r — > oo and, according to Eq. (4.13), v depends 
weakly on d and v oc cj~ 2 (1). The last proportionality shows that the chain velocity strongly 
depends on the interparticle interactions, and v is a decreasing function of r. [Note that the 
permissible values of r are not too close to r±, since our assumptions imply that m&xK <C L 
holds.] 

If K <C max K, then, using the asymptotic formula 40 



F(z + x) 



z x-y 



T(z + y) 

(z — > oo), we reduce the equation for A to the form 



1 

l + 2^(x-y)(x + y-l) 



(4.15) 



^ = 1 + ^ + ^(1-3)- (4-16) 



0! d 6 2 \ d, 

Solving Eq. (4.16) with respect to A and taking into account that in this case 5t <^ 1, we 
derive the formula 

Q x 2 d 



' ~ 1 (d-c)+0 2 c : (417) 
which shows that for 1 <^ K <ti max if the chain velocity does not depend on the chain 
characteristics. Since 1 < A/ R < c/ (d — c), for a given ratio A/R = a (a > 1) the parameter 
c/d must belong to the interval (cr/(cr + 1), 1). On this interval, v is a decreasing function 
of c/d, and according to Eq. (4.17) v — > d\2<yj{G + 1) as c/d — > a /{a + 1), and v — > ^ x as 
c/d -> 1. 



V. CONCLUSIONS 



We have studied the ground state, vibration properties and deterministic transport of 
a chain of charged particles interacting via both the Coulomb interaction and a repulsive 
interaction, which is characterized by the power r. We have shown that a chain with 
equidistant particles has minimum energy for r > r 2 ~ 2.799. For a such chain, we have 
calculated its energy and period as functions of r, derived and analyzed in detail its phonon 
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spectrum, and established that it is stable with respect to small displacements of particles 
from their equilibrium positions only if r > r\ ~ 2.779. 

We have studied also the deterministic transport of a chain in the case when the particles 
are acted upon by both a longitudinal alternating electric field and a force field generated by 
a ratchet potential. We have derived the set of equations that describes the chain transport, 
analyzed the general properties of its solution, and calculated the chain velocity in the 
overdamped limit. If the period of the electric field is large enough, then the chain velocity 
strongly depends on the interparticle interactions. Otherwise, i.e., if the period of the electric 
field is small enough, the chain velocity does not depend on them. 
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APPENDIX A: ANOTHER FORM FOR U (a, (3) 



Using the representation 40 

■ip(z) = lim 



for ip(z) = d\aY{z)/dz (T(z) is the gamma function) and the definition 

OO 1 

C(r,z) = V- - 

of the generalized zeta function £(r, z), we find 



E 

71=1 



+ 



n — a n + a n 



= 2^(1) - if>(a) - i/>(-a) 



and 



E 



i 



n + a\ 



C(r,a) + ((r,l-a). 



(Al) 



(A2) 



(A3) 



(A4) 



Substituting Eqs. (A3) and (A4) into Eq. (2.3) and taking into account that a) = 
— a) + 1/a, 40 we obtain the desired representation 

U (a,P) = -^[ 2 ^(l)-V(«)-V(l-«)] 
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+ ^[2CW + C(r,a) + C(r,l-«)]. 



APPENDIX B: INTEGRAL REPRESENTATION FOR * 3 (#c) AND * 3 (k) 

We proceed from the standard series 39 



^ cos(pz) ( . z 

> = — m 2 sm - 

p=i P \ 2 



(0 < z < 7r). Integrating Eq. (Bl) twice, we find 

l-cos(pz) _ f A r ; , /,., ^ // 

p=l 

and using the definition (3.16), we obtain 



(A5) 



(Bl) 



$ 3 («) = C(3) + J™ dx [ d y ln ( 2 sin I) • ( B3 ) 

To find the integral representation for ^ 3 (k), we use the identity 

g cos(pz) = g cos(2^) + g cos[(2p-l>] (m) 
P =i P P =i 2 P p=i 2p - 1 

and the series (Bl) to obtain 

^ cos[(2p-l)*] 1 / z\ 

g 2 P -l = -2 ln l tan 2j- (B5) 

Integrating Eq. (B5) by the same way as Eq. (Bl) and using the definition (3.17), we find 

7 1 r nK r x ( v\ 

* 3 («) = gC(3) + 2 I dx I ( tan ^) • (B6) 

APPENDIX C: DERIVATION OF EQUATION (4.12) 

Let t k (k = 0, 1, 2K, t = 0, t 2K < t) and t k (k = 0, 1, 2K - 1, t k < t k < t k+1 ) be 
the times at which the positive charges reach the minima and maxima of V(x) during the 
half-period r, i.e., w(t k ) = —A/2 + kd and w(t k ) = —A/2 + kd + c. Solving Eq. (4.11), we 
find 

W (t) = y - (y + f - ^) e -'<"*> (Ci) 
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for t E (t k ,t k ), and 



w(t) = °j - ( 9 j + | - kd - c] e~«) (C2) 



for i G (£fc,£fc+i). At t = Eq. (CI) leads to 



e *fc-fc) = gi/j + A / 2 ~ fC3 ) 
9J8 + A/2 -kd-c' { } 



and Eq. (C2) at £ = t^+i yields 

e *(fc + i-fc) = gg/j + A/2 - fed - c 
9 2 /5 + A/2- kd-d' 

Multiplying Eqs. (C3) and (C4), substituting the result 

5(tfc+1 _ tfe) = 9,/ 5 + A/2 -kd 9 2 /5 + A/2 - kd - c 
9 1 /5 + A/2 - kd - c 9 2 /S + A/2 - kd - d 



(C4) 



(C5) 



into the identity e 5t2K = UI^q 1 e'(* fc + 1- **\ and using the formula nl=i — k) — T(x)/T(x - 
2K) and condition A = 2Kd, we obtain Eq. (4.12). 
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FIG. 1: Part of the equilibrium chain of charged particles. 
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FIG. 2: Plots of Pi(r) and /^(r); for the third solution (3^(r) = oo. The function fair) increases 
from 2 {ri) = ra 4.959 to /3 2 (3) = oo. 
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FIG. 4: Plots of U\{r) and TJ%(r) [figure (a)], and U\(r) and [^(r) [figure (b)]. The functions Us(r) 
and U2(r) are very close to each other [U 3 (r) < U 2 (r) for n < r < 3, f/ 3 (3) = f/ 2 (3) = -l/3\/3, 
and ^(n) — [^3(^1) « 0.001], therefore the plot of Us(r) in figure (b) is not shown. 
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FIG. 5: Dependence of the phonon spectrum Cj{k) on k for r = 4 (curve 1), r = 3 (curve 2), and 
r = ri (curve 3). Since for r = r±, the function |w'(k)| goes to infinity very slowly [as \J\ ln(l — k)\] 
as k — > 1, this feature of the spectrum is not represented on the chosen scale of the k axis. 
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FIG. 6: Illustrative plots of the reduced alternating electric field h(t) [figure (a)], the ratchet 
potential V(x) [figure (b)], and the reduced force field g(x) [figure (c)]. 
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